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FOREWORD 


The  work  reported  herein,  performed  during  the  period  October  1962  — 
January  1965,  was  accomplished  by  the  Cornell  Aeronautical  Laboratory,  Inc, 
(CAL),  Buffalo,  New  York,  for  the  Fluid  Dynamics  Branch,  Seaworthiness 
and  Fluid  Dynamics  Division,  Hydromechanics  Laboratory  of  the  David  Taylor 
Model  Basin  (DTMB),  Washington,  D.  C.  This  research  is  part  of  the  Bureau 
of  Ships  Fundamental  Hydromechanics  Research  Program,  S-R009  01  01, 
administered  by  DTMB  and  performed  under  Office  of  Naval  Research  Contract 
Nonr  3578(00).  The  project  monitor  for  this  program  at  DTMB  was,  initially, 

Mr.  J.  Patton  and,  later,  Mr.  W.  Souders.  Dr.  P.  Crimi  and  Mr.  J.  M.  Grace 
of  CAL  conducted  the  study  under  the  direction  of  Dr.  I.  C.  Statler.  Assistance 
in  digital  computer  programing  was  provided  by  Mr.  H.  Selib. 
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ABSTRACT 


The  flutter  characteristics  of  a  two-dimensional  hydrofoil  near  a 
free  surface  are  analyzed,  with  consideration  given  to  the  cases  of  both 
one  and  two  degrees  of  freedom.  The  flutter  equations  incorporate  lift  and 
moment  as  determined  from  an  exact  linearized  potential  solution.  Free- 
surface  effects  have  been  included  without  approximation. 

Results  indicate  that  for  the  system  having  two  degrees  of  freedom, 
the  free  surface  has  very  little  effect  on  flutter  boundaries,  despite  the 
large  changes  in  hydrodynamic  loading  caused  by  foil-induced  surface  waves. 
When  the  hydrofoil  was  constrained  to  move  only  in  pitch,  however,  flutter 
solutions  were  found  where  none  exist  in  the  absence  of  the  free  surface, 
for  practical  values  of  mass  ratio  and  pitch-axis  position. 
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SUMMARY 


The  study  reported  herein  was  directed  to  analytically  investigating 
the  flutter  characteristics  of  a  two-dimensional  hydrofoil  near  a  free  surface. 
The  problem  was  treated  for  the  hydrofoil  having  either  one  or  two  degrees 
of  freedom.  The  flutter  equations  incorporate  the  lift  and  moment  as 
determined  from  the  exact  linearized  potential  solution  for  the  flow  about 
a  fully  wetted  hydrofoil  oscillating  in  pitch  and/or  plunging.  The  effects  of 
the  free  surface  are  incorporated  in  this  solution  without  approximation. 

Because  the  hydrodynamic  loading  is  a  function  of  two  parameters 
involving  flutter  speed  (i.  e.  ,  Froude  number  and  reduced  frequency),  it  was 
necessary  to  generate  flutter  boundaries  by  a  method  of  cross -plotting. 

This  situation  is  analogous  to  the  one  encountered  in  flutter  calculations  for 
compressible  flow,  in  which  case  loading  depends  on  Mach  number  as  well 
as  reduced  frequency. 

The  hydrodynamic  lift  and  moment  are  greatly  altered  by  the  effects 
of  the  foil -induced  surface  waves  when  the  flow  is  oscillatory.  It  was 
expected,  therefore,  that  a  comparison  of  results  with  those  for  infinite 
depth  would  reveal  large  differences  in  the  flutter  characteristics.  It  was 
found,  however,  that  when  the  system  is  free  to  oscillate  in  both  pitching 
and  plunging,  there  is  very  little  effect  on  flutter  speed  attributable  to  the 
free  surface.  Omission  of  free-surface  effects,  even  for  submergence 
depths  as  small  as  one  semichord,  would  generally  produce  the  correct 
variations  in  flutter  speed  for  a  wide  range  of  parameter  values,  with  the 
results  usually  being  somewhat  conservative. 

The  situation  was  found  to  be  markedly  different  when  the  system 
was  constrained  to  move  only  in  pitch.  The  results  contrast  sharply  with 
those  for  infinite  submergence  depth.  While  in  the  latter  case  flutter  only 
occurs  when  the  pitch  axis  is  forward  of  the  quarter-chord  point,  free- 
surface  effects  cause  flutter,  at  a  relatively  low  forward  speed,  for  the 
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pitch  axis  positioned  aft  of  the  quarter-chord  and  forward  of  midchord.  In 
the  absence  of  the  free  surface,  fluid  density  must  be  much  less  than  that 
of  water  for  flutter  in  a  single  degree  of  freedom  to  occur.  The  results 
obtained  indicate,  though,  that  at  finite  submergence  depths  a  hydrofoil 
would  flutter  in  a  fluid  with  the  density  of  water.  It  is  pointed  out  that  flutter 
of  this  type  could  present  a  problem  when  a  hydrofoil  is  constrained  to  move 
primarily  in  pitch,  as  could  be  the  case  when  it  is  used  as  a  control  surface. 
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1. 


Introduction 


The  similarity  in  the  design  problems  of  aircraft  and  hydrofoil  boats 
has  understandably  led  to  numerous  applications  of  aeronautical  technology 
in  the  analysis  of  hydrofoils.  Of  particular  interest  have  bt  en  the  recent 
investigations  of  hydroelastic  instabilities  (see,  for  example,  References 
1  and  2)  which  correspond  directly  to  flutter  of  aircraft.  The  hydrofoil 
flutter  problem  may  be  contrasted  with  that  of  the  airfoil,  though,  in  that  the 
hydrofoil  is  subjected  to  the  effects  of  a  free  surface.  The  unsteady  motions 
of  a  hydrofoil  cause  surface  waves  wh  ch  markedly  alter  the  hydrodynamic 
loadings. 

Previous  analyses  of  hydrofoil  flutter  have  either  omitted  free-surface 
effects  or  crudely  approximated  them,  because  an  exact  representation  was 
not  available.  The  study  of  hydrofoil  flutter  reported  below,  however, 
incorporates  the  lift  and  moment  as  computed  from  an  exact  linearized 
potential  solution  for  the  flow  due  to  an  oscillating  hydrofoil  near  a  free 
surface.  The  derivation  of  this  solution  is  reported  in  Reference  3. 

This  study  was  directed  to  determining  the  extent  of  the  influence  of 
the  free  surface  on  the  hydroelastic  characteristics  of  a  two-dimensional 
hydrofoil  having  either  one  or  two  degrees  of  freedom  (i.  e.  ,  pitching  and/or 
plunging).  No  attempt  was  made  to  carry  out  a  complete  parametric  study. 
Rather,  selected  system  parameters,  primarily  mass  ratio  and  the  ratio  of 
uncoupled  natural  frequencies,  were  varied,  so  as  to  define  flutter  boundaries 
which  could  be  compared  with  those  that  would  obtain  in  the  absence  of  a 
free  surface. 

In  Section  2,  the  equations  of  motion  are  developed  and  the  method 
used  for  computing  flutter  points  s  outlined.  The  results  of  the  computations 
are  presented  and  discussed  in  Section  3. 
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2. 


Development  of  the  Flutter  Equations 


2.  1  Discussion  of  the  Hydrodynamic  Forces 

As  a  consequence  of  the  study  reported  in  Reference  3,  it  is 
evident  that  the  hydrodynamic  lift  and  moment  acting  on  an  oscillating  hydro¬ 
foil  can  be  significantly  affected  by  the  presence  of  the  free  surface.  The 
extent  of  the  influence  of  the  free  surface  can  best  be  seen  by  reproducing 
certain  of  the  results  of  Reference  3.  Consider  a  hydrofoil  of  infinite  span 
and  chord  2b  oscillating  in  pitching  and  plunging,  as  shown  in  Figure  1. 


FIGURE  1 

Schematic  of  the  System  Analyzed 

The  expressions  for  the  complex  amplitudes  of  the  hydrodynamic  lift  and 
moment  per  unit  span  may  be  written  as  follows  (see  Reference  3),  with 
lift  L  nondimensionalized  by  27T/>U26  and  moment  M  by  npU2^'. 

1  s  aa°)*  T**0*  ikHn(h''a‘ao)  +  H2L  0) 

m*  (’■£  *  t)06*  *  * H**  (,+ir)  a°  <2) 

M  is  the  moment  about  midchord  (positive  nose  up),  /?  a  and  cu  are 

cJb 

defined  through  Figure  1,  and  *s  ^e  reduced  frequency.  The  functions 
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and  //^depend  on  L  .  D  and  the  Froude  number  r  ; 

Each  of  these  functions  reduces,  in  the  limit  of  infinite  submergence  depth, 
to  the  classical  Theodorsen  function.  Thus,  the  deviation  of  any  given 
function  from  the  Theodorsen  result  can  be  considered  as  a  measure  of  the 
effect  of  the  free  surface  on  the  hydrofoil  lift  and  moment.  Since  the 
functional  form  of  each  of  the  functions  is  quite  similar,  it  is  only  necessary 
to  consider  one  example  to  illustrate  the  general  effect  of  the  free  surface. 
As  an  example,  the  real  part  of  Hu  is  plotted  in  Figure  2  versus  reduced 
frequency,  for  a  depth  of  one  semichord  and  a  Froude  number  of  unity.  The 
curve  is  compared  with  the  real  part  of  the  infinite -depth  (Theodorsen) 
function  in  the  figure. 


k 


FIGURE  2 

Comparison  of  H *  for  F*  1  and  D*1  with  the  real  part 
of  the  Theodorsen  function 
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One  can  see  that  the  most  striking  characteristic  of  this  function  is  its 
behavior  near  the  so-called  resonance  condition  (i.  e.  ,  where  fc  =  //V^*). 
Physcially,  the  resonance  point  can  be  considered  as  that  point  where  the 
group  velocity  of  the  wave  caused  by  foil  oscillation  is  equal  to  the  velocity 
of  the  foil.  The  result  of  this  is  an  increase  in  the  amplitude  of  the  waves 
with  an  attendant  significant  change  in  the  velocity  distribution  about  the  foil. 
The  modified  velocity  distribution,  in  turn,  greatly  affects  the  lift  and 
moment. 


The  difference  in  behavior  of  the  hydrodynamic  functions  for 
finite  and  infinite  depth  is  the  motivation  for  this  study,  since  the  marked 
change  of  these  functions  near  the  resonance  point  could  have  an  appreciable 
effect  on  single-  and  two-degree -of-freedom  instabilities.  It  should  be  noted 
that,  in  the  results  presented  below,  flutter  data  corresponding  to  values  of 
reduced  frequency  lower  than  the  resonant  frequency  are  lacking.  Generally, 
there  were  insufficient  lift  and  moment  data  to  adequately  define  a  flutter 
boundary  for  this  range  of  reduced  frequencies.  However,  calculations  for 
selected  values  of  reduced  frequencies  below  resonance  produced  flutter 
speeds  greater  than  those  calculated  for  reduced  frequencies  above  resonance. 

2.  2  Equations  of  Motion 

The  equations  of  motion  for  a  hydrofoil  undergoing  simple 
harmonic  pitching  and/or  plunging  motions  in  the  absence  of  structural 
damping,  when  nondimensionalized,  may  be  written  in  the  form  (see,  for 
example.  Reference  4) 

h0p,  +  ^XCC(X0  -  h0/A,  (^)  ~  -  0  (3) 

hoM'x<*‘ha'»tLr*~  0coM'ra  +~£i(M+2al)  =  O  (4) 

In  the  above  equations,  denotes  the  mass  ratio  f/u  =  m/rr^ob*  )  ,  x&b 

the  distance  of  the  center  of  mass  aft  of  the  elastic  axis,  rab  the  rad  us  of 
gyration  about  the  elastic  axis,  uJf,  the  uncoupled  natural  bending  frequency, 
and  ot)a  the  uncoupled  natural  torsion  frequency.  Combination  of  equations 
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(3)  and  (4)  above  with  (1)  and  (2)  of  Section  2.  1  yields  the  following  form  for 
the  equations  of  motion: 


(m*0 -(<*> h /<*>)*)  -2fl9)  -  *  2iSh]  •O  (5) 

hol(^x'a.+Ai  )  +  ao  [(f*'rJ(1  -(<*>ct/oo)*+  &,)  *  0.  (6) 

The  quantities  f3^  ,  *  ’  Am  and  6t  ,  S2  ,  &3  ,  S¥  are  defined  as: 

_  fH,L  ±  ZOjHf,  I 

A*  -/if 


The  superscripts  (*)  and  (J)  denote  real  and  imaginary  parts,  respectively. 
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If  Equations  (5)  and  (6)  are  to  have  a  nontrivial  solution,  the 
determinant  of  the  coefficients  of  h0  and  oca  must  vanish.  This  requirement 
provides  two  relations  which  must  be  satisfied  at  a  flutter  point,  since  the 
real  and  imaginary  parts  of  the  (complex)  determinant  must  each  equal  zero. 
The  method  used  for  obtaining  solutions  to  these  two  equations  is  discussed 
below. 


2.  3  Method  of  Solution  of  Flutter  Equations 


2.  3.  1  Two-Degree-of-Freedom  System 


The  method  of  solution  used  in  this  report  for  solving 
the  flutter  equations  in  the  case  of  two  degrees  of  freedom  consisted,  first, 
of  specifying  the  system  parameters  D  ,  F  ,  k  , yu,  ,  X>a  ,  and  cl 
and  solving  for  the  values  of  (a^/u>)and  (oo^/u)) that  would  satisfy  the  equations 
in  question.  It  should  be  noted  that  the  specification  of  P  ,  F  ,  and  k.  was 
sufficient  to  determine  the  hydrodynamic  lift  and  moment.  The  numerical 
calculations  were  performed  on  an  IBM  7044  computer.  The  computer 
program  used  is  described  in  Appendix  A  Section  (a). 


The  second  step  used  in  generating  appropriate  flutter 
boundaries  was  to  plot  an  intermediate  graph  of  the  form  shown  in  Figure  3, 
where  the  ordinate  is  given  by  *  Fk.  .  It  is  necessary  to  obtain, 

first,  this  intermediate  plot  rather  than  the  usual  plot  of  vs‘ct)h/^ 

from  the  solutions  U^fu),  cJ^/u) ,  and  the  input  data,  since  the  latter  plot 
would  consist  of  a  family  of  curves  with  Froude  number  as  the  parameter. 
These  curves  could  not  be  regarded  as  flutter  boundaries,  since  both  the 
ordinate  and  the  parameter,  Froude  number,  are  functions  of  velocity.  The 
purpose  of  using  the  intermediate  plot  with  ordinate  is  to  remove  this 

difficulty.  The  situation  is  analogous  to  the  one  encountere  d  in  the  flutter 
analysis  of  an  airfoil  in  compressible  flow,  where  loading  is  a  function  of 
reduced  frequency  and  Mach  number.  The  final  form  of  flutter  boundary  is 
generated  from  Figure  3  by  constructing  a  cross -plot  of  (JjbuJa  vs. 
for  constant  values  of  where  the  ordinate  is  given  by 

u  _  r 


bu>n 


ou_ 
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Flutter  boundaries  with  /c  as  an  independent  variable 

were  also  constructed.  These  were  generated  by  using  an  intermediate 

graph  similar  to  Figure  3,  except  that,  for  this  case,  Froude  number  is 

held  constant  and  is  the  parameter  of  the  family  of  curves.  Data  from 

a  series  of  these  graphs,  each  for  a  constant  Froude  number,  were  then 

U 

obtained  for  a  given  value  of  thereby  allowing  a  plot  of  vs.  /ul 

to  be  generated. 


2.  3.  2  Single -Degree-of-Freedom  System 

The  flutter  equations  for  a  single-degree-of-freedom 
system  can  be  deduced  from  Equations  (5)  and  (6).  For  example,  a  foil 
which  has  infinite  stiffness  in  torsion  has  the  following  flutter  equations: 

+ -Zf  =  O  (7a) 

=  O-  (7b) 

For  the  case  in  which  there  is  infinite  stiffness  in  the  bending  mode,  Equation 

(6)  gives 

/+rj [t  -(rter/e*)9]  +  &i  *0  (8a) 

6Z  -  O'  (8b) 

Since  only  flutter  in  torsion  was  treated,  comments  concerning  the  method  of 
solution  refer  to  Equation  (8).  However,  the  methods  of  solution  of  Equations 

(7)  and  (8)  are  analogous.  The  computer  program  used  to  solve  Equation  (8a) 
and  (8b)  are  described  in  Appendix  A  Section  (b). 


As  noted  previously,  the  specification  of  the  parameters 
D  ,  F  ,  and  k  are  sufficient  to  determine  the  hydrodynamic  functions. 
Therefore,  Equation  (8b)  may  be  solved  for  a  once  0  ,  f  ,  and  Jc  are 
specified.  Equation  (8a)  may  then  be  solved  for  (u^/cujusing  the  value  of  CL 
calculated  from  (8b)  and  the  specified  values  of  .  The  solution  for 

can  be  written  as 


(9) 
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A  plot  of  vs.  a,  can  then  be  generated,  where  the  ordinate 

expressible,  from  Equation  (9),  in  the  form 


f!  •  rtl 


1  + 


Si  V/i 


(10) 


^a\ 

It  was  found  that  a  more  complete  definition  of  flutter 
boundaries  could  be  obtained  with  the  hydrodynamic  data  available  if  , 

rather  than  a,  ,  were  made  the  independent  variable.  Thus,  intermediate 
plots  of  vs.  O/  were  first  formed,  for  specified  values  of  D  and  F  , 

with  as  parameter.  As  an  example,  one  of  these  intermediate  plots 

is  shown  in  Figure  4  for  D*f  and  F  */  .  Next,  cross -plots  of  vs. 

,  for  a  specified  value  of  cu  ,  were  constructed,  with  Froude  number 
as  parameter.  Finally,  values  of  were  selected  and  the  curves  of 

U/bu)a  vs.  generated  from  the  latter  plots. 
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Figure  4  EXAMPLE  OF  PLOTS  USED 
FOR  THE  SYSTEM  HAVING 


TO  GENERATE  FLUTTER  BOUNDARIES 
ONT  DEGREE  OF  FREEDOM 
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3. 


Discussion  of  Results 


3.  1  Flutter  in  Two  Degrees  of  Freedom 

The  computations  with  both  pitching  and  plunging  degrees  of 
freedom  included  were  directed  primarily  to  obtaining  an  indication  of  the 
magnitude  of  free-surface  effects,  rather  than  to  generating  a  complete 
parametric  study.  Thus,  submergence  depth  was  maintained  at  one  semi¬ 
chord,  while  a  single  value  of  0.  2  for  mass  imbalance  Xa  was  used  and 
r*  was  kept  at  a  value  of  0.  35. 


In  Figures  5,  6  and  7,  dimensionless  flutter  speed  (S/ba^is 
shown  plotted  against  for  values  of  Ou  of  -0.  1,  -0.  2  and  -0.  3, 

respectively.  On  each  figure,  the  flutter  boundaries  for  of  0.  75  and 

1.  35  are  shown,  together  with  the  curve  for  infinite  submergence  depth. 
The  value  of  is  1.3  for  all  these  boundaries. 


Most  surprising  about  these  results  is  the  similarity  between 
the  curves  for  finite  and  infinite  submergence  depths.  Figure  5  shows  no 
appreciable  effect  of  the  free  surface  on  flutter  speed.  Figures  6  and  7 
indicate  that,  if  the  pitch  axis  is  moved  forward,  the  free  surface  has  some 
stabilizing  effect,  but  the  shift  in  the  curves  can  hardly  be  termed  substantial. 
It  might  be  reasoned  that  the  surface  waves  which  are  generated  would  tend 
to  stabilize  the  system  by  carrying  off  energy,  but  the  effects  would  be 
expected  to  be  much  larger  because  the  hydrodynamic  loading  is  so  greatly 
affected  by  the  free  surface. 


In  Figures  8,  9  and  10  are  shown  flutter  boundaries,  obtained 
by  cross-plotting,  with  p,  as  abscissa,  for  of  0,  2,  0.  5  and  0.  75, 

respectively.  Again,  each  figure  shows  the  curves  for  u of  0.  75  and 
1.  35  together  with  the  one  for  infinite  depth.  Pitch -axis  position  Qj  is  -0.  2 
for  all  three  of  these  figures.  The  free  surface  is  again  seen  to  have  very 
little  effect  on  the  flutter  speed.  As  in  the  previous  plots  with  0^/u^.as 
independent  variable,  the  surface  waves  have  some  stabilizing  influence, 
particularly  for  =  1.  35  and  for  the  two  smallest  values  of  “V«i- 
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Figure  7  Ujba)A  vs.  fuJ*  for  <x>  -  -0.3 


BB  - 1 629-S-3 


Figure  10  U/buJ^v s.  /<,  for  =  0.75 


It  would  appear,  then,  that  despite  the  large  effect  of  the  free 
surface  on  the  hydrodynamic  loading,  the  flutter  of  hydrofoils  having  two 
degrees  of  freedom  can  be  adequately  predicted  without  taking  account  of 
the  free  surface  for  submergence  depths  as  small  as  one  semichord. 
Omission  of  free-surface  effects  would  generally  provide  conservative 
results,  and  the  over-all  effects  of  various  system  parameters  would  be 
correctly  obtained. 

3.  2  Flutter  in  a  Single  Degree  of  Freedom 

In  contrast  with  the  situation  for  two  degrees  of  freedom,  the 
dominance  of  free-surface  effects  in  the  hydrodynamic  loading  was  found  to 
greatly  alter  the  flutter  characteristics  of  a  hydrofoil  with  just  one  degree 
of  freedom.  A  system  of  the  type  being  analyzed  here,  in  the  absence  of 
free-surface  effects,  is  extremely  resistent  to  flutter  when  one  of  the  two 
degrees  of  freedom  is  removed.  A  two-dimensional  hydrofoil  at  infinite 
submergence  depth  will  not  flutter  at  all  if  motion  is  allowed  only  in  plunging, 
and  if  it  is  just  free  to  pitch,  flutter  only  occurs  when  the  pitch  axis  is  ahead 
of  the  quaj  ter-chord  point  (i.  e.  ,  a.<  y-  )  and  the  mass  ratio  ju.  is  much 
greater  than  would  be  encountered  by  a  hydrofoil  (see,  for  example,  Refer¬ 
ence  5).  The  addition  of  free-surface  effects  was  found  to  alter  this  situation 
markedly. 


Consideration  was  limited  to  flutter  in  pitch.  It  appears  that 
flutter  in  plunging  is  possible,  but  it  would  occur  over  a  very  restricted 
range  of  reduced  frequencies  for  which  lift  and  moment  data  were  not 
available.  The  results  of  the  computations  are  shown  in  Figures  11  and  12, 
where  the  dimensionless  flutter  speed  is  plotted  against  for 

selected  values  of  pitch-axis  position  and  .  The  values  for  a,  selected, 

i.  e.  ,  -.112,  -.  201  and  -.  266,  were  not  chosen  completely  arbitrarily; 
judicious  choices  for  the  values  of  cb  greatly  facilitated  the  definition  of 
flutter  boundaries  because  of  limitations  in  the  amount  of  hydrodynamic 
data  available.  Again,  computations  were  only  performed  for  a  depth  of 
one  semichord.  It  appears,  however,  that  flutter  characteristics  similar 
to  those  obtained  would  be  evident  at  other  finite  depths. 
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The  results  of  Figures  11  and  12  indicate  that  the  effects  of  a 
free  surface  are  extremely  destabilizing  for  flutter  in  pitch  alone.  At  finite 
submergence  depth,  flutter  occurs  when  the  pitch  axis  is  aft,  rather  than 
forward,  of  the  quarter -chord  point.  Instead  of  there  being  a  minimum  in 
/^^a,  below  which  no  flutter  occurs,  as  is  the  case  for  the  system 

becomes  progressively  less  stable  as  decreases.  Generally,  the 

reduced  frequencies  for  flutter  were  found  to  lie  quite  near  the  critical 
values.  This  unusual  flutter  instability  may  therefore  be  attributed  to  the 
singular  behavior  of  the  hydrodynamic  loading  near  the  resonance  condition. 

Instabilities  of  this  type  could  conceivably  be  significant  in 
some  hydrofoil  applications.  For  example,  a  hydrofoil  used  as  a  control 
surface  could  have  a  pitch  frequency,  as  determined  by  the  stiffness  of  the 
control  linkages,  which  is  much  less  than  its  bending  frequency.  The  foil 
would  tnen  effectively  have  only  a  pitching  degree  of  freedom  and,  so,  would 
be  subject  to  the  type  of  flutter  predicted  here. 
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APPENDIX  A 


DESCRIPTION  OF  THE  COMPUTER  PROGRAM 

(A)  "Hydrofoil  Flutter  Program  2  Degrees  of  Freedom" 

Written  by  H.  Selib  at  the  Cornell  Aeronautical  Laboratory  ,  Inc. 

The  purpose  of  the  IBM  7044  Fortran  IV  program  generated  was 
to  obtain  a  solution  to  the  equations  of  motion  (3)  and  (4)  of  Section  2.2. 
These  equations  were  solved  to  obtain  the  frequency  ratio  (  )  for 

given  values  of  the  input  parameters  namely,  D,  k,  F,  a,  u.  ,  r*  »  z  # 
3a,  '  3h  •  anc*  *he  hydrodynamic  functions  H^,  H^y  ^2M’  aS 

defined  in  Reference  3. 

As  noted  above  the  program  was  written  in  Fortran  IV  and  consists 
of  a  main  program  of  length  ^  2330  (number  of  octal  locations). 

The  inputs  are  as  follows: 


CARD 

FORMAT 

1 

16, F8,  p 

NC  (number  of  cases  for  one  set  of 

D,  F,  k,  a  values);  a 

2 

9F8.  5 

D,  F,  k 

3 

9F8.  5 

rh1l,  ih1l;Rh2L,  ih2L;RH1m, 

ihim;  rh2M;  ih2M 

4 

9F8.  5 

up  to  Card  NC  +  4:/^  ra  »  Set 

The  program  print  out  lor  the  case  of  two  degrees  of  freedom  flutter  consists 
/•■/ A  uJl 

of  values  of  — —  ,  ■ -  and  products  of  these  quantities  with  F  and  k  which  are 

CO  OU 

of  interest.  In  addition  quantitites  Lj,  Ly  Ly  Ly  M^f  My  My  M^  are 
included  in  the  print  out.  These  quantities  can  be  related  to  quantities  defined 
on  page  5  of  this  report  as  follows: 

23  BB-1629-S-3 


I 


I 


L,  a 

III 

L  2  ~  ' 

uz~  Az 

L3  =-lS* 

L  4=-2<S^ 

M  *-Sz 

The  program  listing  is  contained  on  the  following  pages. 
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SlttFTC  ZZFLTR  LIST.RLF 

C  HYDROFOIL  FLUTTER  PROGRAM  -  TWO  DEGREES  UF  FREEDOM 

C  REVISED  -  5/20/63 

DIMENSION  AAMUI 1 00  )  . XXA ( 100  I , RRA ( TOO  I , GGA 1 10 0 ) , GGH (l 00 ) 

CUMMUN  A AMU, XX A, RK A, GGA, GGH 

COMMON  D,F,AK,A,NC,R1L,XI1L,R1M,XI 1M , R 2L , X 1 2L , R2M , X I 2M , XL l , XL2 , 

1  XL3,XL4,XMl,XM2»XM3, XM4, AMU »XA,RA,GA,GH,V1 ,V2 ,AKl,AK2,AK3, 

2  AK4,Cl,C2,C3,C4,B0,Hl,B2,0ISC,LL,SQD,SIGl,SiG2,SlHl,SlH2, 

3  WA12,MA22, WH12,WH22,MAM1,MAW2,WHRA1,WHWA2, VBWAl , VbMA2,KHHAl 

4  tRH8A2,CHBAl,CH»A2,FBMAVl,FBtaAV2 
LuoO  FORMAT!  16, F8. 5) 

l-ul  FORMAT  ( 9F8 • 5  I 
1^02  FORMAT!////  ) 

Ijj3  FORMAT!  1H1,42X33HTW0  DIMENSIONAL  HYDROFOIL  FLUTItR  //  ) 

i^o4  FORMAT l 33X3HD  *F6.3,10X7HA  »F  3 . 3 ,  1 0X4HXA  *F6.3  /33X3HF  =F6.3, 

1  IJX7HMU  =F8 • 3, I 0X4HGA  *F6.3  /33X3HK  *F6. 3 , 1 0X7HRA*»2  * 

2  F8 • 3 , 10X4HGH  *F6. 3  //  ) 

1  0  u  5  FURMAT  (41X  38 HD IbCKlMINANf  IS  NEGATIVE  -  NO  SOLUTION  III  ) 

1  -  J6  FORMAT (21X12H(mA/W)  l*  *2  *F12.3,  5X 12H{ WA/W) 2**2  -F12.3,  2X1H-,  2X 
1  11HNO  SOLUTION  III  ) 

1-07  F ORMAT !21X12H(WA/W)2»*2  *F12.3,  5X48H ( WA / w )1 • *2  IS  POSITIVE  -  NO  S 
1ULUT1UN  FOR  (WA/W)2  //  ) 

i jG8  FORMAT!  33X12H|Wm/W) 1»«2  *F l 2 . 3 , 5X2 3HNU  SOLUTION  FOR  (WH/W)1  //  I 
10  09  FORMAT  (23XllH(WA/W)  1  *F12.3  /20X 1 1H ( WH/MA ) 1  =F1?.3  /20X 

1  1 1 H ( V/ HHA ) i  *F12«3/20X11H(FBWA/V) 1  *F 1 2 . 3/2GX 1 1  HR! H 3/BA ) 1  * 

2  F12.3  /20X11HI IHQ/BAI1  *F12.3  ) 

UID  FORMAT!  33X12HI*m/W)2*»2  *F 1 2 . 3 , 5X2 3HNO  SOLUTION  FOR  tWH/w)2  //  I 
1  1  1  FORMAI  12JXUHIWA/M)  1  »F12.3,30X11H(  WA/W)2  *F12.3  /2DX 

1  11H|WH/WA)1  =F i2. 3,  30XliH(WH/WA)2  *F 12. 3/2DXI 1H( V/HWAI 1  * 

2  F12.3, 1GX11H! V/BWAJ2  *F 12 . 3/20X UH ! FBtfA /V ) 1  *Fi2.1,3^X 

3  1 1 H ! FBh A/ V  )  2  =F12.3/20X11HR!H0/BA) l  *F 1 2 . 3 , 3GX 1 1 HR  I  HO/BA ) 2  * 

4  Fi2.3/20X11HI!H0/!1A)1  *F  12 . 3,  30X 1 1M I  ( HO/BA  )  2  *F12.3  I 

1-12  FORMAT  !  48X24HAEROOYNAM IC  COEFFICIENTS  // 30X4HL 1  =FB.3,5X4HL2  *FH.3 

1  ♦ 5X4HL 3  =F8 • 3, 5X4HL 4  *F8.3  /30X4HM1  *F8.3,5X4HM2  =F8.3, 

2  5X4HM3  *F 8 • 3 , 5X4HM4  *F8.3  //  ) 

2  READ  100 J, NC, A 
READ  10b 1 , 0, F • AK 

3  READ  lOOi.RlL.Xl l L , R2L , X I 2L , R 1 M, XI IM,R2M,XI2M 

4  XL  1  ^  1.0*2.J»XI 1L/AK 
XL2  *  -2.0*R1L/AR 

X L 3  *  -A*!  XI2L-2.0*A*XI1U/AK-2.0*R2L/AK**2 
XL4  =  -I 1.0^R2L-2.0*A*RILI/AK-2.G*XI2L/AK*»2 
6  XMl  =  -A-!  XI  IM*2.0*A*XI1L  )/AK 
XM2  =  !RiM*2.3*A*RlU/AK 

X  M3  =  .125*A*«2>! A*  (  X I 1M-XI2LJ-.5*! XI 2M-4,«A**2*XI IL ) ) / AK* ( R2M^2 • • 
1  A*R2L)/AK*#2 

6  XM4  =  !A*!1.3*R2L  -R 1M ) - . 5* < l . 0-R2M ) -2 . J* A*»2* Rl L ) / AK* ! X I 2M« 2. 0* A* 

1  XI2LI/AK#*2 

DO  7  1*1, NC 

READ  100 l, AAMUI 1 ) ,XXA( I),RRA< I ) , GGA ( I ) , GGH! I ) 

7  CONTINUE 

DO  40  1*1, NC 
AMU  *  AAMU (  I ) 

X A*XXA ( I ) 

8  RA  *  RRA!  I  ) 

GA  *  GGA!  I  I 
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GH  *  GGH( I  ) 

9  VI  -  GH» AMU-XL2 
V2  =  GA*AMU*RA-XM4 

10  AK1  *  V2-GH*XM3 

AK2  *  XL1»V2^XM**VI*AMU»XA»(XM2+XL4)*XL3*XM2*XM1*XL4 
AK3  *  GA+GH 
AK4  *  GA*XL l-Vl 

11  Cl  *  -(XLUGA (VI) 

C2  =  -XLl*XM3^Vl*V2*AMU*XA*Umj*XA*XL3^XMl)*XL3*XMl-XM2*XL4 
C3  a  1 .0-GA*GH 
C4  *  XM3*GH*V2 

12  BO  *  C  l * AK  3  +  C 3* AK4 

B 1  =  C2*AK3«-C4»AK4-Cl*AKl-C3*AK2 
B2  =  “(C2*AKl*C4*AK2J 
IFIB0H20,  13*13 
120  BO  =  -BO 
HI  *  -B I 
B2  *  “B2 


13 

DISC  *  Bl**2-4.j 

*B0*B2 

IF(LL) 14, 15. 14 

14 

wRirt 

GO  TU  16 

(6 

*  1002 

) 

15 

MK 1  T  L 

(  6 

,  1003 

) 

16 

MR  11  L 

(6 

*  1004)0, A, XA,F,AMU,GA,AK,KA,GH 

1  60 

mR  1 1 1 

(6 

• 1012)XLl,XL2,XL3,XL4,XMlf XM2 , XM3 • XM4 

17 

LL  *  LL  + 1 

IF(LL“3) 19, 18, 10 

18 

LL  *  0 

19 

1  F  (  01  SC  )  20  *  2  1  *  2 1 

20 

MR  1 1 1 

GO  f Li  4C 

(  6 

,  1005 

) 

21 

SUL)  *  SUKT  (DISC) 

S 1 G 1  *  ( -Bl-SQD) / ( 2 

•  0*BQ ) 

SIG2  *  (-Bl*SQDI/(2 

•0*B0) 

S  1  Hi  *  (S1GI»C1+C2)/(SIG1*C3+C4) 

SIH2  *  (SIG2*CI*C2)/ISIG2»C3*C4> 

22  MA12  *  l.O-SIGl/ ( AMU*RA) 

MA22  *  l.0-SIG2/( AMU*RA) 

NH12  *  l.O-SIHl/AMU 
mH22  *  1.0-SIH2/AMU 

23  I  F  (  M  A 1  2)24*24,25 

24  MKIJL  (6  , 1006)MA12,MA22 

GO  I U  40 

25  NAMl  >  SORT  ( M A 12  ) 

IF(MA22>26,26,260 

26  MRirt  (6  *1007) MA22 

260  IF(MH12)27,28,28 

27  MRITE  (6  *  1008 ) MHl 2 

MHMA l  *  99999.0 

GO  TO  29 

28  MHMAl  =  SORT  (MH12)/MAWI 

29  VBMAl  =  1.0/ ( AK*WAW1  ) 

FBMAV1  =  F#AK*MAMl 

KH8A1  *  (AMU*XA*XMl)*(AMU»RA»(MAl2-l.0)-XM3)fXM2»(MAl2«GA“XM4) 
CHBA1  *  I  AMU*XA^XH1)*(WA12*GA-XM4)“XM2»( AMU*RA«( MA12-1. ))-XM3) 

30  I F ( MA22 ) 31  *  3 1  *  32 
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31  irfRITk  (6  •  1009)  WAMl  ,  WHWA  l ,  Vtt WA  l  »  FB WA  V 1  »  KHBA 1  ,  LHBA l 

SO  TO  40 

32  MAW2  *  SORT  ( W A22 ) 

I F ( WH22  > 33.34,34 

33  *RIU  (6  » 1010)  MH22 

MHWA2  *  99999 • 3 

SO  TU  3b 

34  mHmA2  *  SORT  (MH22I/MAM2 

35  VBWA2  *  1 .0/ ( AK» WAW2 l 
FBWAV2  -  F  *AK *WAW2 

RHBA2  *  (AMU*XA*XM1)*(AMU*RA*(WA22-1.0)-XM3)*XM2*(WA22»GA-XM4) 
CHBA2  *  { AMU*XA*XM1 )*( WA22*GA-XM4>-XM2*f AMU*RA* C WA22-1 . J ) -XM3 > 

36  NRlTb  (6  ,  1011)  WAWi  ,  WAW2,  WHNA  i  ,WHMA2,  VB  WA  1 ,  VliW  A2  , 

1  FBMAV1 1  F-8WAV2«  RHBA l,RHBA2tCHBAlfCHBA2 

40  SONT I NUfc 
SO  TU  2 
tNU 


27 


BB-1629-S-3 


(B)  "Hydrofoil  Flutter  Program  ...  1  Degree  of  Freedom  - 

Pitching" 

Written  by  H.  Selib  at  the  Cornell  Aeronautical  Laboratory,  Inc. 

The  purpose  of  this  IBM  7044  Fortran  IV  program  is  to  solve 
Equations  (8a)  and  (8b)  of  Section  2.  3.  2  for  (  ^5  and  a,)  once  the  values  of 

i  00 

D,  F,  k,  have  been  prescribed.  The  specification  of  D,  F,  k  is 

sufficient  to  determine  the  hydrodynamic  functions  H^,  H^,  and  H^, 
^e^nec*  Reference  3,  and  involved  in  Equations  (8a)  and  (8b). 

In  the  following  the  main  program  of  length  g  1013  (number  of 
octal  locations)  was  written  in  Fortran  IV.  The  subroutine  clear  of  lengthg 
27,  used  to  zero  common  storage,  was  not  written  in  Fortran  IV  but  is  a 
"MAP"  program. 

The  inputs  are  as  follows: 


CARD 

FORMAT 

1 

9F8.  5 

2 

M'rct 

2 

9F8.  5 

D,  F,  K.  If  D  =  99999999 

a  new  value  for 

start  with 

3 

9F8.  5 

RH1L’  IH1L:  RH2L'  IH2L' 
RH2M*  IH2M* 

RH1M’  ihim* 

The  program  print  out  for  the  single  degree  of  freedom  system  consists  of 

values  of  a,  ,  oolJ A,  &  — —  as  functions  of  D,  F,  k  and  Lur*  • 

^  *  iS 

The  program  listing  follows  on  the  next  page. 
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HBFTC  ZZ1FLT  LIST, REF 

C  HYDROFOIL  FLUTTER  PROGRAM  -  UNE  DEGREE  OF  FREEDOM 

COMMON  AMURA2,NLIN,D,F,AK, F IL , G 1L , F 2L , G2L , F 1M , G IM , F 2M, G2M , B, C , 

1  DISC. IFLU, JFLG.KFLG, Ai,A2,WAl2,MA22,WdGl,MBG2,UBMl ,UBM2, 

2  to AL . NA2 
EQUIVALENCE  (D.ND) 

1  READ  1000 « AMUKA2 
IjOO  FORMAT  (9G8.5) 

NLiN  =  0 

2  CALL  CLEAR(D,MA2l 
READ  1000.0.F.AK 
IF1NU.EQ. 999999991  GO  TO  1 

READ  1000.FlL.GlL.F2L.G2LtFlM.GIM.F2M.G2M 

3  H  =  (  1.u*F2L-F1M*2.0*C2L/AK)/I4.0*FIL) 

C  =  (F2M*2.0*G2M/AK-1.0)/(4.0*FIL > 

DISC  =  B**2*C 

IF(DISC.GT.O.J)  GO  TO  4 
1 F LG  *  2 
GO  IU  9 

4  I  FLO  =  1 

A 1  =  B  +  SQR  T ( 0 1  SC ) 

A2  *  B-SQKT (DISC) 

MA12  *  1.0  + (A  1**2* ( 1.0*2.0»GIL/AK)4AI*HGIM-G2L)/AK*2.0*F2L/AK«*2> 
l  ♦.125-.5«G2M/AK*F2M/AK»«2)/AMURA2 

mA22  *  l.GMA2**2*<  1 .0^2 . 0*GIL/ AK )  *A2*  I  (  G1M-G2L)  /AK*2.0*F2L/ AK*«2 ) 
l  ♦.125-.5*G2M/AK+F2M/AK»*2)/AMURA2 

5  IMWA12.GT.0.0)  GO  TO  6 
JFLG  »  2 

toAl  =  -SQRTI-WA12) 

GO  TO  7 

6  JFLG  =  1 

toAl  =  SORT l WA12 ) 
toBGl  *  to  A 1  *  AK  *F 
JQW1  *  l.O/(AK*WAl) 

7  IF(toA22.GT .0.0)  GO  TU  8 
KFLG  =  2 

to A2  =  -SORT ( -WA22 I 
GO  TU  9 

8  KFLG  *  l 

mA2  =  SQK1 (WA22) 
toHG2  =  WA2*AK*F 
JBm2  =  1 .0/ ( AK*m A2 ) 

9  I F ( NLIN.EQ.J)  WK I TE ( 6, 1001 ) AMURA2 

1JUI  FORMAT  t IH1 , 52X26HHYDRUFU I L  FLUTTER  IN  PITCH  //53XI4HMU*RALPHA**2  * 

1  L12.5  //14X1HD, 10X1HF, iOXIHK, 14X1HA, 14X4HtoA/to,  8X 

2  12HWA»SQKT(B/G>,  7X6HU/B*toA, l JX6HB»»2*C  ) 

GO  TU  (11*  1  j  )  . I FLG 

10  toRlTE(6, 1002ID.F. AK.01SC 

1,02  FORMAT  IF  17.2, F 13.5, F 10. 5.F8 1.5) 

ML  IN  *  NL I N+ 1 
GO  Tu  17 

11  GO  TO  (  12.  13) .JFLG 

12  MRITEI6, 1003)0. F.AK.Al. MAI, WBG1.UBW1.DI SC 

1 j03  FORMAT (Fl7.2,FlJ.b,Fl0.5,Fl5.5,F16.b,F18.S,Fl5.5,F17.5) 

NLIN  -  NLINM 
GO  T  (J  14 
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13  riRi  T  fc  (  6f  1004)D#F.AK,A1,WA1,  DISC 

1304  FORMAT  (F17*2»F13«5»F10*5»F15«5»F16*5|F53*5) 
NLIN  *  NLINM 

14  CU  TO  ( 1 5» 16) « KFLG 

15  mRI T  t ( 6 *  1005 )  A2 f  MA2 • MBG2 • UBM2 

L j05  F0RMAT(F55.5,F16.5» F 10. 5# FI 5. 5) 

NLIN  *  NLIN+l 
CO  TU  17 

16  *RlTfc<6, 1005)  A2*WA2 
NLIN  *  NLIN+I 

17  IFINLIN.GE.5J)NLIN  *  0 
CO  TO  2 

bNU 
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